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The lens surface profile is derived based on the instantaneous focal length versus the lens radius 
data. The lens design based on instantaneous focal length versus the lens radius data has many 
useful applications in software assisted image focusing technology. 
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I. INTRODUCTION 

The software assisted image focusing is an emerging 
technology which is expected to replace the traditional 
methods of autofocusing in image manipulation devices, 
such as digital cameras and mobile phones. Unlike the 
traditional methods, where the focusing of images is 
done by mechanically movable parts, the software as- 
sisted technology produces focused images by processing 
it through specialized image reconstruction algorithm. 
The transition from mechanical to software assisted im- 
age focusing can be attributed to the (1) demand for 
thinner and lighter products by customers, and (2) the 
advancements in manufacturing process for faster and 
more power efficient digital signal processors. With aut- 
ofocusing by mechanically movable parts, the demand for 
thinner and lighter products is becoming a top hurdle for 
manufacturing process. On the other hand, the advance- 
ments in more power efficient and faster digital signal 
processors make software assisted image focusing tech- 
nology ideal for satisfying customer's demand for thin- 
ner and lighter image developing products such as digital 
cameras and mobile phones to name a few. 

At the heart of software assisted image autofocusing 
technology is the specialized image reconstruction algo- 
rithm permanently coded into the built in digital sig- 
nal processor. The actual layout of the code base for 
image reconstruction algorithm varies among different 
manufacturers and many manufacturers do not disclose 
their algorithms to public as they constitute a trade se- 
cret. The image reconstruction algorithm can be codified 
based on instantaneous focal length versus the lens ra- 
dius data(author?) Once this specialized image re- 
construction algorithm is adopted for the system, a lens 
must be designed so that its output matches the instan- 
taneous focal length versus the lens radius data, which 
information was assumed and used as input to the image 
reconstruction code base. 

In this work, a formula for the lens surface profile is 
presented. The derivation of lens surface profile is solely 
based on the instantaneous focal length versus the lens 
radius data; and therefore, the result is expected to find 
useful applications in software assisted image focusing 



technologies. 



II. INSTANTANEOUS FOCAL LENGTH DATA 



Instant focal distance p as function of lens radius y 
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Figure 1: Instantaneous focal length versus the lens radius, 
where both are measured in meters. The normal incidence is 
assumed for the incoming light rays. 

Alexander and Lukyanov have recently filed for a 
patent which deals with image reconstruction algorithm 
with applications in software assisted image focusing 
technology. In their proposal, they claim to have ob- 
tained an optimal image processing solution, which is 
expected to be a significant improvement over the prede- 
cessor(author?) [2]. Behind their optimization is the in- 
stantaneous focal length versus the lens radius data illus- 
trated in Fig. [U which result assumes a normal incidence 
for the incident light waves. In the figure, f3 denotes the 
instantaneous focal length and y is the lens radius. Each 
of the twelve segmented curves can be represented by the 
quadratic polynomial 



Pi (y) = a^y"^ + hy + a, 
with coefficients (ai, bi, Ci)given by 

ai = -313.07, bi = 0.0235, ci = 0.0035137034, 



(1) 



a2 = 534.53, 62 = -0.2472, C2 = 0.003527877626, 
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03 = -309.02, 63 = 0.0818, C3 = 0.0035088062, 
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a4 = 536.05, 64 = -0.3275, C4 = 0.0035493232, 



as = -306.12, 65 = 0.1182, cg = 0.003502912672, 



as = 539.03, be = -0.3891, cg = 0.003569538239, 



07 = -303.68, 67 = 0.1463, C7 = 0.003496845208, 



as = 542.21, bs = -0.4417, cg = 0.003589312176, 



ag = -301.27, 69 = 0.1695, cg = 0.00349080193, 



aio = 545.96, 610 = -0.4895, cio = 0.003609151039, 



an = -298.81, 611 = 0.1893, cn = 0.003484870978, 



ai2 179.08, 612 = -0.0474, = 0.003469596542, 

where the subscript i denotes the ith curved segment in 
Fig. [TJ The curve fitting was done by Hnear regression. 
The physical lens, whose output satisfies the instanta- 
neous focal length versus the lens radius data defined in 
Fig. [11 is one of the variants of lens with extended depth 
offocus(author?) ^MM- With Eq. (H]) , I shall solve for 
the lens surface profile whose output matches the instan- 
taneous focal length versus the lens radius data defined 
in Fig. [H 



III. THE LENS SURFACE EQUATION 
A. Derivation 

When a ray of light passes across media of different 
refractive indices, its path is governed by the Snell's law, 



sm ( 



ng sm ( 



(2) 



as illustrated in Fig. [21 Here, = {cj) and ne = 
ng (uj) are frequency dependent refractive indices with uj 
denoting the angular frequency of the light. The param- 
eters (j) and 6 represent the angle of incidence and angle 
of refraction, respectively. 

If N denotes the normal vector to the local point y = 7 
on the curve x = h{y) , then it can be shown 



||-Nx(-ei)|H||-N||||-ei||sin 
and the expression for sin becomes 
INxeil 



A^sinc 



sm = 



N 



N^\\N\ 



(3) 



Similarly, the expression for sinO may be obtained by 
considering vectors A, B, and C of Fig. [51 The vectors 
A, B, and C satisfy the relation, 

A + B = C. (4) 

In explicit form, vectors A and B are defined as 

A = -762, B = (/3-a)ei, (5) 

where 62 is the unit basis for the y axis. With Eqs. ([4]) 
and 1(5]), the vector C becomes 

C 61 -762. (6) 

The vector cross product N x C is given by 

N X C = (/? - a) N X 61 - 7N X 62 

and its magnitude becomes 

||N X C|| = ||(/3- a)N X ei - 7N x 62II = NC sinO, 

(7) 

where N = ||N|| and C = ||C|| . Utilizing Eq. ©, C may 
be expressed as 



1/2 



and the Eq. ([7]) is solved for sin 6 to yield 

||(/3 - a) N x 61 - 7N X 62! 



snif 



N 



(/3-«f+7' 



1/2 



(8) 



Refracted ray 



x=h(y) 



(y=0, x=0) 
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Figure 2: Illustration of Snell's law. 

Insertion of Eqs. ^ and ^ into the Snell's law of Eq. 
^ gives 



where ei is the unit basis for the x axis. 
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ng 



|(/3 - a) N X 61 - 7N X 62I 
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(9) 
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By definition, the normal vector N satisfies the relation, 

where g [x, y) is a function whose gradient gives N, 

„ do do dh 
N = V.g = —ei + —62 = ei - -^^2. 

ox oy oy 

Because N is the normal vector at the location 
{x — a,y — ^) ,\ write 



N = ei - — 

oy 



62- 



(10) 



The following vector cross products are valid. 



N X ei = ei X ei — — — 
dy 

N X 62 = ei X 62 - 

oy 



62 X ei, 



J/=7 



62 X 62, 



J/=7 



With these replacements, Eq. I|T2|) gets re-expressed in 
form as 



dx 
dy 



y 



P-x-'-^\{P-xf + y^ 



1/2- 



(13) 



How is the instantaneous focal function, /3, restricted? 
The (} in Eq. lfT3|) is restricted so that the expression for 
dx/dy does not blow up. Equation l|13p is well behaved 
if and only if the denominator satisfies the condition. 



(3-x 



ng 



(/3 - xY + y' 



1/2 



^0. 



To solve for /?, I shall first rearrange the previous expres- 
sion to get 



P-x^^ 

710 

Squaring both sides. 



(/3 - xf + y^ 



1/2 



(/3 - x) V 4 (/3 - + 4?/'' 
ni ni 



where Eq. (|TOl) was used to replace N. Since ei x 6i 
62 X 62 = 0, the previous relations reduce to 



N X 61 = — 
oy 



63, NX62=63, 



(11) 



V=7 



where 63 is the unit basis for the z axis of which satisfies 
the relation, 

61 X 62 = 63, 62 X 61 = —63. 

Insertion of Eq. (fTTjl into Eq. ^ gives 



and regrouping the terms, I find 

f 



The resulting expression can be solved for /? to yield 



I3^x± 



n,py 



(14) 



Equation l|T4l) defines the restriction for the instanta- 
neous focal function. 



ne dh 



y=i 



■ 7 



y=7 



(f3-af+r 



1/2 ' 



which expression can be rearranged to yield 

dh 
dy 



1/2 ' 



(12) 



where a and 7 are constants of which are depicted in Fig. 

m 

For Alexander and Lukyanov's optical element, the in- 
stantaneous focal function (3 = f3 (y) in Eq. lfT2|) is as 
defined in Fig. [TJ The 7 for the y axis is not anything 
special, of course. Any y belonging to the domain of h 
satisfies the Eq. I|12p . The generalization of Eq. (fT2|) for 
all y belonging to the domain of h is done by making the 
following replacements: 



dh 
dy 



V=7 



dh dx 
dy dy' 



B. Lens surface profile 

The profile of axially symmetric lens about its optical 
axis is obtained by solving the initial-value differential 
equation, Eq. (flSl) . 



dx 
dy 



y 



where x {yo) = xq is the initial condition to be specified 
and the instantaneous focal function /? satisfies the con- 
strain defined in Eq. (fT4|) . Without loss of generality, one 
may choose x {y = yo = 0) = for the initial condition 
and the lens profile satisfies the differential equation. 



dx 
dy 



iP,-xf + y^ 



1/2 ' 



(15) 



x (0) = 0, (3i^x± 
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where the index i in Pi comes from the fact that the input 
specification defined in Fig. [T] is piece wise continuous 
over range of x. The domain for each /3i is given by 

/3i : 0<y< 0.00019182692, 
(32 ■■ 0.00019519231 <y < 0.00027259615, 
/33 : 0.00027596154 <y< 0.00033317308, 
Pi : 0.00033653846 <y< 0.00038701923, 
P5 ■■ 0.00039038462 <y< 0.00043413462, 

Pe : 0.0004375 <y< 0.00047451923, 
Pr : 0.00047788462 <y< 0.00051153846, 
Ps : 0.00051490385 <y< 0.00054855769, 
Pg : 0.00055192308 <y< 0.00058221154, 

Pio : 0.00058557692 <y< 0.0006125, 
Pn : 0.00061586538 <y< 0.00064278846, 

Pi2 : 0.00064615385 <y< 0.00067307692. 

The differential equation ifTSj) has been solved using the 
Runge-Kutta method(author?) p]. The Runge-Kutta 
routine has been coded in FORTRAN 90 and the result 
for the case where — 1 and ng = 1.5311 is provided in 
Fig. [31 The physical lens may be obtained by revolving 
the curve about the x axis. Since the light ray is directed 
in the positive x direction, as illustrated in Fig. IH it 
implies that the image sensor should be embedded inside 
the lens for the case where ne > n^. 
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Figure 4: Location of the image sensor. 



Reversing the values for two refractive indices, i.e., 
(n0 = 1.5311 and ne = 1), the lens surface profile be- 
comes as illustrated in Fig. [5l Again, the physical lens 
may be obtained by revolving the curve about the x axis. 
Since the light ray is directed in the positive x direction, 
the case < represents the situation where light is 
exiting the lens medium. For this configuration, where 
ne < n^, the image sensor should be placed external to 
the lens medium, as illustrated in Fig. [H 



Cross-sectional view of lens 



Cross-sectional view of lens 
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Figure 3: The lens cross-section, ng > n^. 



Figure 5: The lens cross-section, ng < n^. 
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Image sensor 
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Detailed cross-sectional view of lens boundary 
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Figure 6: Location of the image sensor. 



3.5e-()6 
3e-()6 
2.5e-()6 
2e-()6 
^ 1.5e-()6 
le-()6 
5c-()7 


-5e-07 

0.0001 0.0002 0.0003 0.0004 0.0005 0.0006 0.0007 
Lens radius v, [m] 

Figure 7: Plot of a; — a;i for the case where = 1 and 
ne = 1.5311. 
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The plots of lens surface cross-section, Figs. [3] and [5l 
deceptively portray as if lens surface profile is represented 
by parabolic class of curves. To show graphically that 
this is not the case, the numerical data solutions obtained 
via Runge-Kutta method to graph Figs. [3] and [5] were 
linearly regressed to obtain 



xi = aiif" + hiy^ 
+ fiy + 9i, 



ciy 



diy^ 



eiy 



(16) 



X2 = a2y^ + b2y^ + C2y^ + d2y^ + e2y^ 

+ /2y + 52, (17) 

where xi is the polynomial curve fit for Fig. [31 X2 is the 
polynomial curve fit for Fig. O and the coefficients are 
given by 

ai = -8 X 10l^ bi = 2x 10", ci = -1 X 10^ 
di = 52523, ei = 403.1, /i = 3 x lO^^, 
51 = -3 X 10"^ 



02 = 5 X 10", 62 = -1 X 10" 
d2 = -29916, 62 = -264, /2 
<?2 = 2 X 10-9 



, C2 = 1 X 108, 
= -2 X 10-4, 



If xi represents a perfectly fitting polynomial functions 
for the curve plotted in Fig. O then one should expect the 
difference x—xi is a constant, where x is the plotted curve 
in Fig. [31 Similarly, if X2 represents a perfectly fitting 
polynomial functions for the curve plotted in Fig. [H then 
one expects the difference x — X2 is a constant, assuming 
X now is the plotted curve in Fig. [5l Contrarily, if x — xi 
(or X — X2) is not a constant, then the polynomial Xi (or 
X2) cannot be a perfectly fitting polynomial function for 
the curve plotted in Fig. [31 (or Fig. [5]). 



Detailed cross-sectional view of lens boundary 
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Figure 8: Plot of a; — a;2 for the case where = 1.5311 and 
ne = 1. 



The results are shown in Figs. [71 and [8l respectively for 
the cases where ng > and ne < n^. It clearly shows 
that X — xi or X — X2 are far from being constants. This 
indicates that the surface cross-sectional profile of lens 
is not a simple parabolic curve as deceptively portrayed 
Fig. [31 (or Fig. [5]). Instead, the magnification of the 
surface reveals series of kinked segments which must be 
attributed to the discrete continuous curve segments in 
instantaneous focal length (/3) versus the lens radius (y) 
data shown in Fig. [H 



IV. CONCLUDING REMARKS 

At the heart of software assisted image focusing tech- 
nology is the specialized image reconstruction algorithm, 
which is permanently coded into the built in digital sig- 
nal processor. The algorithm is often codified basing on 
the instantaneous focal length versus the lens radius data 
as the initial input. The software assisted image focusing 
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system therefore requires a specially designed lens whose 
output generates the instantaneous focal length versus 
the lens radius data. In this work, a formula for the lens 
surface profile has been presented. The derived lens for- 
mula generates a unique surface profile for the lens based 
on the instantaneous focal length versus the lens radius 
data. The lens design based on instantaneous focal length 
versus the lens radius data makes this result well suited 
for software assisted image focusing technology. 
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